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1 INTRODUCTION 



ABSTRACT 

Extending previous studies, we derive generic predictions for lower order cumulants and their 
correlators for individual tomographic bins as well as between two different bins. We derive the 
corresponding one- and two-point joint probability distribution function for the tomographic con- 
vergence maps from different bins as a function of angular smoothing scale. The modelling of 
weak lensing statistics is obtained by adopting a detailed prescription for the underlying density 
contrast. In this paper we concentrate on the convergence field k and use top-hat filter; though 
the techniques presented can readily be extended to model the PDF of shear components or to in- 
clude other windows such as the compensated filter. The functional form for the underlying PDF 
and bias is modelled in terms of the non-linear or the quasilinear form depending on the smooth- 
ing angular scale. Results from other semi-analytical models e.g. the lognormal distribution are 
also presented. Introducing a reduced convergence for individual bins, we are able to show that 
the tomographic PDFs and bias for each bin sample the same functional form of the underlying 
PDF of density contrast but with varying variance. The joint probability distribution of the conver- 
gence maps that correspond to two different tomographic bins can be constructed from individual 
tomographic PDF and bias. We study their dependence on cosmological parameters for source dis- 
tributions corresponding to the realistic surveys such as LSST and DES. We briefly outline how 
photometric redshift information can be incorporated in our computation of cumulants, cumulant 
correlators and the PDFs. Various approximate results for cumulants and their correlators are pre- 
sented. Connection of our results to the full 3D calculations is elucidated. Analytical results for 
inclusion of realistic noise and finite survey size are presented in detail. 

Key words: : Cosmology- Weak-Lensing - Methods: analytical, statistical, numerical 



Following the first weak lensing measurements dBeacon. Refregier&EUisI l2000l : | Wittm an et a il I2000I: Kaiser, Wilson & Luppind [2000 ; 



IWaerbeke et all20 00) , the field of weak lensing has witnessed a tremendous progress in all fronts (see lMunshi et al J d2008h for a review). Currently, 
in terms of cosmological observations, weak lensing plays a role complementary to both Cosmic Microwave Background (CMB) studies and stud- 
ies involving large scale structure (LSS) surveys. The ability of weak gravitational lensing to reveal cosmological information, particularly the 
dark energy equation of state is considerably enhanced by the inclusion of tomographic information. The impotance of weak lensing has spurred 
tremendous progress on the technical front in terms of specification and control of systematics. There are many ongoing and future weak lensing 
surveys such as the CFHlQ legacy survey, the Pan-STARRS^ and the Dark Energy survejQ and further in the future, the Large Synoptic Survey 



1 http://www.cfht.hawai.edu/Sciences/CFHLS/ 

2 http://pan-starrs.ifa.hawai.edu/ 

3 https://www.darkenergysurvey.org/ 
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Telescopfl Joint Dark Energy Mission or JDE1V0 that will map the dark matter and dark energy distribution of the entire sky in unprecedented 
details. In particular, owing to the large fraction of the sky coverage and tighter control on systematics as well as dense sampling of source galaxy 
populations it will be soon possible to study gravity induced non-Gaussianity with extreme accuracy. The gravity induced non-Gaussianity is 
typically probed using re al space correlation functions as well as in the harmonic domain using their harmonic counterparts i.e. the multispectra 
(see e.g. |Pen et al] d2003l) ). These correlation functions provide a set of tools to go beyond the usual power spectrum analysis. The higher-order 
correlation functions are important not only to break the parameter degeneracy inherent in power spectrum analysis (e.g. between the amplitude of 
the matter power spectrum as and the matter density paramet er Qm) but also to understand error-es timates of lower-order correlations functions. 
Starting with the study of the three-point correlation function dvillumserj[l996l ;l Jain & Selia j 1997b higher order statistics of weak lensing shear, 
convergence or flexions are now well understood from a theoretical point of view. 

The power spectrum of density perturbations remains the most commonly used statistic in many cosmological studies. Weak lens- 
ing surveys probe the non-linear regime and are sensitive to non-Gaussianity which can not be probed using only the two-point cor- 
relation function or its harmonic analo g the power spectrum. The statistics of shea r or convergence probe the statistics of underlying 
mass distribution in an unbiased way dJain, Seliak & White] |2000|: IMunshi & JairJlioOll ; iMunshil I2OO0I ; iMunshi & Jair] |2000| : IValageasI 2000; 
IValageas, Munshi & Barbeil 1200 5; Takada & Whitell2003l ; iTakada & Jainll2004l) . sensitive to nonlinear evolution due to gravitational clu stering . 
Various analytica l schemes from perturbati v e calc ul ations to halo models have been employed to model t h e weak lensing statistics iFrvl dl984 ); 
Schaeffej dl984l); iBernardeau & SchaeffeJ dl992l); ISzapudi & Szalayl dl993l Il997ri; IMunshi et all dl999h : IMunshi. Coles & Melottl Jl999albh : 
Munshi, Melott & Colesl Jl999h : IMunshi & Colesf d200q 12002. l2003h : ICoorav & Sethi J2002h ). In addition to studying the statistics in projec- 



tion they have also been studied in 3D using phot ometric redshifts. This approach can further tighten the constraints on e.g. the neutrino 
mass as well as the dark energy equation of state (Heavens 2003; Heavens, Refregier & Hevmans 200~(il; llfeavens, Kitching & TavloJ [20061 ; 
iHeavens, Kitching & Verde] |200 71; ICastro et all 120051 ; [Kitching et al. 2008). Tomographic techniques have also been employed as an intermedi- 
ate strategy between projected surveys and 3 D mapping (Hu 1999; Takada & Jain 2004, 2003; Sembo loni et a ll2008h. 

In this paper we extend previous results jjain. Seliak & White 2000; Muns hi & JanT.20 01 ; Munshi 2000; M unshi & Jainl200d ; IVaiageasl2000h 
on projected surveys by analysing the entire one-point PDF and the two-point PDF with tomographic information.The PDF contains information 
about the correlation hierarchy to an arbitrary order; the correlation hierarchy of the convergence field is directly related to that of the underlying 
mass distribution. We employ a generating function formalism that relies on hierarchical ansatz on smaller angular smoothing scales and on 
perturbative results on larger scales. We define a reduced convergence for each bin and show that the different bins sample the same underlying 
PDF and bias functions (to be defined later) for the density contrast. The entire joint two-point PDFs for different pairs of redshift bins and 
individual PDF for each bins can be constructed from the PDF and the bias associated with individual bins because the joint PDF is factorisable in 
terms of the individual PDFs, bias and cross-correlations among various bins and different angular scales. We will show that individual redshift- 
resolved tomographic maps can be used to map out the PDF of the underlying mass distribution for a wide range of variance. This underlying PDF 
of the density contrast can be used to recover the tomographic PDF with the use of just two individual variables Kmin and the reduced variance for 
each bin; both of these variables are uniquely determined by the geometry and matter content of the Universe. The results are applicable not only 
to the PDFs as determined under hierarchical ansatz but also for other well motivated approximations for PDF such as the lognormal distribution. 

Recent cosmologi cal observations favour an accelerating Univ erse. This implies existence of energy of unknown nature (dark energy) which 
has negative pressure dAmendola & TsuiikawdPioiol ; IWangi 2010). Current data continues to be consistent with dark energy being a non-zero 
cosmological constant. Though many other alternative dark energy candidates have been consider which are consistent with data as well, e.g. 
quinessence, k-essence, spintessence. Different dark energy models can be classified according to the equation of state of of the dark energy 
component «7x- For quintessence model dwx/dz > while for k-essence models dwx/dz < 0. There are many complimentary probes for dark 
energy, the distance-redshift relation of cosmological standard candles; Cosmic Microwave Background Anisotropy; volume redshift relations 
using galaxy counts; the evolution of galaxy clustering; weak lensing, etc. The different methods to probe dark energy are complementary to 
other and can provide important consistency check. Weak lensing surveys are particularly suitable for dark energy studies. All major weak lensing 
surveys has dark energy as their one of prime science driver. We will use the techniques developed in this paper to study two different dark energy 
model and compare the predictions against those of standard ACDM model. The methods presented here are complementary to the usual Fisher 
matrix based approach that rely on two-point correlation functions or the power spectrum as it includes non-Gaussian information upto order . 

This paper is organised as follows. In §2 we introduce our notation and present some standard results. In §3 we link the lower order statistics of 
weak lensing convergence to that of the underlying density distribution. In §4 we briefly review the hierarchical ansatz in the context of generating 
function formalism. In §6 we discuss the lognormal model in the context of weak lensing statistics. In §7 we derive the PDF and bias for various 
tomographic bins. The results are quite generic and can be used for arbitrary source redshift distribution. Finally the §8 is left for discussion of our 
results. In an appendix we outline how in the context of tomographic binning the evolution topological estimators such as Minkowski Functionals 
can be studied using the lognormal distribution. 



4 http://www.lsst.org/llst_home.shtml 

5 http://jdem.gsfc.nasa.gov/ 
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2 NOTATION 

The statistics of the weak lensing convergence /t(ft) represents that of the projected density contrast 5(x) along the line of sight. In our analysis 
we will consider a small patch of the sky where we can use the plane parallel approximation or small angle approximation to replace the spherical 
harmonics by Fourier modes. The 3-dimensional density contrast S along the line of sight when projected onto the sky with the weight function 
ws(r, r s ) gives the weak lensing convergence in a direction Cl which we have denoted by k(Ci): 

Single Source Plane : K (ft) = f dr ws(r) *(r, ft); uj s (r,r a ) = |^ Amu" 1 M^If ~ r) ; Kf in (r s ) = - P ws(r, r s )dr. 

Jo z c dA(r s ) Jo 

(1) 

Here dA(r) is the angular diameter distance at a comoving distance r. The subscript s in ws(r, r s ) refers to a single source plane. We have also 
introduced a parameter « mln which will be useful in parametrization the PDF and represents the minimum value of the convergence k\ Ho is the 
Hubble parameter and a represents the scale factor. The comoving radial distance is denoted by r. For a distribution of sources represented by 
p s {z) we can write the projected convergence «(ft) as follows: 

u>s{r,r s ) = ^-^-n M a' 1 (r)-^-dA(r) f dr s p s (z) ^- -^—r- ; p s {z) = exp(-z/z ). (2) 

Z c Tig J r ar s u,A\r s ) ZZq 

In a tomographic analysis the source population is divided into several redshift bins and each of which is treated separately. The contribution 
from the individual bins are taken into account when computing the cumulants or the cumulants correlators. It is also possible to compute the 
cross-covariance of these redshift bins. The convergence (ft) from i-th tomographic bin can be expressed as: 

Tomography: K W (ft) = f w (i) (r)S[r, ft]; W(i){r) = ^^2-fi M ^-a _1 (r) d A (r) f dr s p s (z)^- dA ^ s ^ (3) 

Jo Z C Ui Jmax{r, ri } ar s «A(,fsJ 

The "bin average" of the source population is denoted by n, and is defined accordingly fii = J r,+1 dr 3 p s (z)dz/dr s . We will consider different 
bin sizes and source distributions. To incorporate the photometric redshift error we can write 



rr i+ i f 

Photometric Redshift Errors : (r) = r / dr 

J n 



Zh) 



F K (r',r). (4) 



Here Fk = [SK(r — r')/5 , K(f)5 , K(f')] with Sk(r) = sinh(r), r, sin(r) for open, flat and closed geometries. The probability distribution p h 
signifies posterior probability distribution of redshift given a photometric redshift of zu- In our calculation we will need to define a new variable 
K mm (or K™j n for tomographic bins) which will be useful later: 

KT n (r s )= drw s (r,r s ); K™ n = / dr w (i) (r); (5) 
Jo Jo 

In evaluation of k™" we use the following approximate form for the window to 4 (r): 



3 H 1 
W(i)(r) w Ars-—§-QM—dA{r)p 3 (z(ri)) 

Z C Tli 



dz 



dr s 



d A {rr - r) 



Using these results it is easy to see that k™" = K™ m (ri). 

We will adopt two example survey configurations to make definitive calculation. For DES we will take zo = 0.3 and for LSST we will take 
20 = 0.4. The range of source distribution that we consider for each survey is z s = 0.2 — 1.6. The bin-size we take is Az s = 0.2. The constant 
fig is set by imposing the normalized condition J 00 dz p s (z) = 1. For our purpose we have n g = 1.2 x 10 7 ng (n g specifies the galaxy number 
density per square arc-minutes). We will vary n g from few galaxies per arcmin 2 to tens of galaxies per arcmin 2 . The noise power spectrum Cf in 
terms of the intrinsic ellipticity -yf = 0.1 is expressed as Cf = jf /n g . 

Next we consider the lower order cumulants for individual bins as well as projected catalogs. These results will be eventually be useful for the 
construction of the entire PDF and bias. The particular cosmology that we will adopt for numerical study are specified by the following parameter 
values: fi A = 0.741, h = 0.72, fi 6 = 0.044, S! C dm = 0.215, n M = tt b + ficDM,n s = 0.964, w Q = -l,w a = 0,cr 8 = 0.803, Q„ = 0. 

In addition to the ordinary ACDM model we will also use two dark energy models in our study. The angular diameter distance for a dark 
energy dominated comsology with dark energy equation of state Six can written as: 

d A (z) = cHo 1 ^ dz' [n M (l + z f + n K (l + z'f + tt x f(z)]- 1/2 (7) 
Jo 

Here fix denotes the dark energy component and Ok = 1 — — fix. The function f(z) parametrizes the time-dependence of the dark energy 
density and f(z = 0) = 1. For dark energy with constant equation of state t»x = px/px we have f(z) = (1 + z) A( - 1+ulx \ The ACDM is a 
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limiti ng case when wx = — 1 or f(z) = 1. In general for an arbitrary dark energy equation of state can be represented as l lWang &"G anrnavich 
1200 ll) Wx(z) = 1/3(1 + z)f'(z)/f(z) — 1. The popular parametrization is given by u)x(z) = wo + Wiz/(1 + z). We will consider two different 
models: (i) constant equation of state wo = —0.95; and (ii) with evolving equ ation of state inx(z) = — 1 + 2/(1 + z). 

For the computation of the power spectrum we use the scaling ansatz of jPeacock & Doddslll994) . The ansatz consists of postulating a non- 
local mapping A-Kk z P(k) — f n i [Atv kfPi (fej)] of linear power spectrum Pi(ki) at a wavenumber ki to nonlinear power spectrum P(k) to another 
wave number k. The wave numbers k and ki are related by an implicit relation ki = (1 + Auk 3 P(k))~ 1 ^k. The functional form for f„i is 
determined from numerical simulations (see Sj6]for more related discussions). 

The evolution of the linear power spectrum in a dark energy dominated model can be charetcerised using a function g(z) i.e. Pi(k, z) = 
[g(z) /(l + z)] 2 Pi(k, z = 0). Where g(z) can be expressed as: 

5 A 00 1 -4- z' 

g (z) = -n M {i + z)E( z ) j dz' r^A]3 : = >/n M (i + + n K (i + «) 2 + W (8) 

We will use these expressions to compute the variance of smoothed convergence field k(6o) as a function of source redshift and smoothing radius 
6q. We will use top-hat smoothing window WthO^o) f° r ° ur study. 

In Figure (TJ we plot the parameter K m i n as a function of redshift for different cosmologies (left panel). We also show the number distribution 
of source galaxies (right panel). 



3 LOWER ORDER STATISTICS OF TOMOGRAPHIC CONVERGENCE MAPS 

Using these definitions we can compute the projected two-point correlation function in terms of the dark matter power spectrum Pg (k, r) (Peebles 
1980, Kaiser 1992): 



/ 



d A (r) 



WlniWo). (9) 



Here 812 is the angular separation projected onto the surface of the sky (cos \6i2\ = Hi ■ ^2) and we have also introduced 1 = dA(j")kj_ to 
denote the scaled projected wave vector; uii(r) are the weak lensing projection weights for the ith photom etric bins defi ned in Eq.l[3]l; to include 
photometric redshift errors we simply need to use u>i{r) defined in Eq.©- Using Limber's approximation dLimberiri954l) the variance of ku)(0q) 
smoothed using a top-hat window Wth(#o) with a radius Go can be written as: 

The variance is plotted for different redshift bins in Figure l[2}. Similarly the higher order moments of the smoothed temperature field relate 
(n p (9o)) to the 3-dimensional multi-spectra of the underlying dark matter distribution B v (Hui 1999, Munshi & Coles 2000): 



/ 4 . r s , #H)W r d 2 h r d 2 \ 2 r d% w n . 



(2tt 



\d A (r)' cU(r)' d A (r)' ^aW 'l2h=o 

The subscripts li — represent the delta function foQZ U)< We will use these results to show that it is possible to compute the complete 
probability distributi on function of Kg) from the underlying dark matter probability distribution function. Details of the analytical results presented 
here can be found in dMunshi & ColesfcOOOh . A similar analysis for the higher order cumulant correlators of the smoothed convergence field relating 
(fli)/c?U (&2)}c with multi-spectra of underlying dark matter distribution B p+q can be expressed as dSzapudi & Szalavll99"7l ; lMunshi & Colesl 

booolbooa) : 

xexp(ie 12 -l 3 )Bi(-A T ,-4 T ,--4 T ,r) . (13) 

\d A (r) d A (r) d A {r) /£i;=o 

We will use and extend these results in this paper to show that it is possible to compute the whole bias function b(> k), i.e. the bias associated 
with those spots in convergence map which k is ab ove certain threshold (which acts as a generating function for these cumulant correlators) from 
the statistics of underlying over-dense dark objects dMunshi. Coles & M elon* 1999^13). 
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4 HIERARCHICAL ANSATZE 

The spatial length scales corresponding to small angles are in the highly non-linear regime of gravitational clustering. Assuming a "tree" model 
for the matter correlation hierarchy in the highly non-linear regime, one can write the general form of the iVth order correlation function Q as 
(Peebles 1980, Bernardeau & Schaeffer 1992, Szapudi & Szalay 1993): 

£< 3) (n, ra, r 3 ) = Q s ($° (n, r 2 )£< 2) (n, r 3 ) + cyc.perm.); (14) 
e< 4) ( ri ,--- ,r 4 ) = 7? a (d 2) (ri,r2)ef (ri,r 3 )^ 2) (n, r 4 ) + cyc.perm.) + i? b (Cl 2) (ri, r 2 )e?(r 2 , r 3 )Cf (r 3 , r 4 ) + cyc.perm.). (15) 

In general for correlation functions of arbitrary order are constructed by taking a sum over all possible topologies with respective amplitudes 
parameters Qn.ci, which in general will be different: 

(p-i) 

&\ ri ,...r p ) = E II d 2) (r*, ri ). (16) 

a, p — trees labcilings edges 

(2) ~ 

To simplify the notation we will use Q (ri, r 2 ) = £i 2 and £ 2 for its volume average over a volume v. It is interesting to note that a similar 
hierarchy develops in the quasi-linear regime in the limit of vanishing variance (Bemardeau 1992); however the hierarchical amplitudes Q p , a 
become shape dependent functions in the quasilinear regime. In the highly nonlinear regime there are some indications that these functions 
become independent of shape, as suggested by studies of the lowest order parameter Q 3 = Q using high resolution numerical simulations 
(Sccocimarro et al. 1998). In Fourier space such an ansatz means that the hierarchy of multi-spectra can be written as sums of products of the 
matter power-spectrum: 

Bf^k^ka.ka)^;^^ = Q 3 (P 5 (k 1 )P 5 (k 2 ) + P 5 (k 2 )P 5 (k 3 ) + P 5 (k 3 )P 5 (k 1 )); (17) 
B< 4) (ki,k2,k 3 ,k 4 ) Efci= o = Ra P 4 (ki)P 4 (ki +k 2 )Pi(ki + k 2 +k 3 ) + cyc.perm. + R b ( P <5 (ki)P 4 (k2)Pi(k 3 ) + cyc.perm. (18) 

In general for p-the order poly-spectra (ki , . . . k p ) we can write: 

(p-i) 

B< p) (k 1 ,...k p )= Y, Qp.° E II p *( k » k ;)- ( 19 ) 

a, p — trees labcilings edges 

Different hierarchical models differ in the way they predict the amplitudes of different tree topologies. Bernardeau & Schaeffer (1992) considered 
the case where amplitudes in general are factorisable, at each order one has a new "star" amplitude and higher order "snake" and "hybrid" 
amplitudes can be constructed from lower order "star" amplitudes (see Munshi, Melott & Coles 1999a,b,c for a detailed description). In models 
proposed by Szapudi & Szalay (1993) it was assumed that all hierarchical amplitudes of any given order are degenerate. Galaxy surveys have been 
used to study these ansatze. Our goal here is to show that weak-lensing surveys can also provide valuable information in this direction, in addition 
to constraining the matter power-spectra and background geometry of the universe. We will use the model proposed by Bernardeau & Schaeffer 
(1992) and its generalization to the quasi-linear regime by Bernardeau (1992, 1994) to construct the PDF of the weak lensing field Km . We express 
the one-point cumulants as: 

<4)> c = (3g 3 )c««] = s«(4))? (4 ) > c = (i2R a + 4P b )ci l) [4l = <4)>« m 

where we have introduced the following notation: 

C;\ q [[Je {r)f +q - 2 \Je 12 {r)]] = P ^S^^[Jo (r)] p+q ~ 2 [Jo 12 (r)]dr; (21) 

Jo <Ia ( r ) 

WoW] - / ^ (^)) WkGft,). [Je X2 {r)] = / ^Ps (^) W^o)e X p(l. 6 12 ). (22) 

The normalised cumulants for convergence in the i-th bin are denoted by (skewness) S3 1 ' and (kurtosis) and are plotted in Figure (3)- Eq.l|20t 
was derived by Hui (1998) in the context of weak lensing surveys. He showed that his result agrees well with the ray tracing simulations of Jain, 
Seljak and White (1998). Later studies extended this result to the entire family of two-point statistics such as cumulant correlators (Munshi & 
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0.5 1 1.5 0.5 1 1.5 



Figure 1. The parameter K m i n is plotted as a function of redshift z s in the left panel for various background cosmologies. The right panel shows the source density 
distribution (not normalised) for the two different surveys. The lines along the x-axis denotes the positions of the tomographic bins considered in our analysis. Notice 
that the parameter fc m i n do not depend on smoothing angular scales and only depend on the depth of the survey as well as on redshift distribution of source population. 
We consider two different dark energy models along with ACDM cosmology. The curve qCDM correspond to wo = —0.95 and the model wOwl correspond to an 
evolving dark energy model w(x) = —1 + z/(l + z). 

Coles 1999, Munshi & Jain 1999). 

(4 ) (A 1 )« w (fi 2 )) c = 2Q 3 Cf ) [J eo (r)J 9l2 (r)] = ^ (23) 

<K?o(fii)Ky)(fia))c = (3i?a + 6R b )Ci^[J eo (r)Jg 12 (r)] = C 3 \Cf > [Jg (rf Jg 12 (r)] = cg 3 1 ) ) (4 ) )^^ (l) (n 1 )^(fi 2 )) c , (24) 

(4 ) (n 1 )4 ) (n 2 )> c = 4JfcCl y) [^(r) 3 Ji ia ^ (25) 

(4 ) (f2i)K U) (n 2 )) c = (24S a + 36S b +4S c )C<i 3) [Je Q {rfje 12 {r)] = C^C^lJe^rf Je 12 {r)] = <4)>^ W (&)><£26) 

(4 ) (n 1 )4 ) (f2 2 )) c = (12S a +6S b )d; f) [Je (rfjo 12 (r)] = C^J eo (rf^ (27) 

where Cp q denotes the cum ulant correlato rs for the underlying mass distribution. These results essentially employ the small angle approximation 
or Limber's approximation dLimbeJ 1954b that are routinely used in computation of higher order cumulants in many cosmological contexts. Other 
approximations such as the Born approximation that we use have been verified by testing against simulations. 

In a related, but slightly different context, these lower order statistics can also be helpful in probing the pressure bias as a function of scal e 
for the study of thermal Sunyaev-Zel'dovich (tSZ) effect or its cross-correlation against tomographic weak lensing maps jMunshietaT]|2011bl) . 
The thermal Sunyaev-Zel'dovich effect probes the line-of-s ight integral of electronic p ressure fluctuations. Cross-correlating (frequency-cleaned) 
y maps from ongoing CMB experiments such as PlanckH Planc k Collaborationf 2006) against the weak lensing tomographic maps can provide a 
redshift resolved picture of reionization history of the Universe. The signal-to-noise will however decrease with increasing order of these statistics. 
Indeed the study of PDF or bias that we undertake next will essentially combine information from all orders. 

In Figure (O and Figure ^ we have plotted the variance and the lower order S p parameters respectively as a function of smoothing scale 8q. 



5 THE GENERATING FUNCTION 

In a scaling analysis of the count probability distribution function (CPDF) the void probability distribution function (VPF) plays a fundamental 
role. It can be related to the generating function of the cumulants or S p parameters, 4>(y) (White 1979, Balian & Schaeffer 1989) : 

P v (0) = exp(-N<r(N c )) = exp f - . (28) 

Where P v (0) is the probability of having no "particles" in a cell of of volume v, N is the average occupancy of these "cells", and N c = All 
statistical quantities correspond to underlying density contrast S. The VPF P v (0) is a special case of the count probability distribution function 
(CPDF) P V (N). The VPF P v (0) contains information of the entire CPDF P V (N) 

The VPF is meaningful only for a discrete distribution of particles and can't be defined for smooth density fields such as S or n(9o). However 

6 http://www.rssd.esa. int/index.php?project=SP 
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Figure 2. The plots shows the variance in convergence {k?{9q)) c as a function of smoothing angular scales 9q. A top-hat window has been assumed. The curves 
from top to bottom correspond to various tomographic bins. The redshift bins correspond to Az s = 0.2 and covers a range z s = 0.2 — 1.4. The curve qCDM 
correspond to wq = —0.95. The higher curves correspond to the deeper redshift bins. See text for more details. 



the scaling functions er(y) and cj>{y), defined above in Eg. J281l. u(y) = —<t>(y)/y, are very useful even for continuous distributions where they 
can be used as a generating function of one-point cumulants or S p parameters: <j>(y) = Sp/p\y p . The function <f>(y) satisfies the constraint 

s s L =si = 1 necessary for proper normalization of PDF. The other generating function which plays a very important role in such analysis is the 
generating function for the vertex amplitudes v n associated with nodes appearing in the tree representation of higher order correlation hierarchy 
(Q3 = v 2, R a = v\ and Ri, — v- A ). In practice it is possible to work with a perturbative expansion of the vertex generating function G(t). In 
terms of the vertices this is defined as: Q(t) = 53^=o( — l) n ^n/«-!- However in the highly nonlinear regime a closed form is used. A more specific 

model for Q(t), which is useful to make more specific predictions (Bernardeau & Schaeffer 1979) is given by G(t) = ^1 + t/ko^J . We will 
relate K a with other parameters of scaling models. While the definition of VPF does not involve any specific form of hierarchical ansatz it is to 
realize that writing the tree amplitudes in terms of the weights associated with nodes is only possible when one assumes a factorisable model of 
the tree hierarchy (Bernardeau & Schaeffer 1992) and other possibilities which do not violate the tree models are indeed possible too (Bernardeau 
& Schaeffer 1999). The generating functions for tree nodes can be related to the VPF by solving a pair of implicit equations (Balian & Schaeffer 
1989), 

The above description has been limited to the level of constructing one-point PDF. A more detailed analysis is needed to include the effect of 
correlation between two or more correlated volume element which will provide information about bias and cumulant correlators. The bias b(8) 
can be introduced the following expression for the joint or two-point PDF: 

P (S 1 ,5 2 )d5 1 dS 2 = p(5i)p(5 3 )(l + fo(<5i)£i26(<52))d<Sid<52 (30) 

The function r(y) - so metimes denoted by /3(y) in the literature - plays t he rol e of a generating function for the factorized cumulant correlators 
c pi ( C pq = c pi c qi) feemardeau & Schaeffd[T9^ : lBemardJl99l[T994l) : r(y) = E^Li V P Cpi/P l - We wil1 next consider two different 
regimes; the quasilinear regime valid at large angular scales and the highly nonlinear regime valid at smaller angular scales. 



5.1 The Highly Non-linear Regime 

bb The PDF p(S) and bias b(S) can be related to their generating functions VPF cj>(y) and r(y) respectively by following equations (Balian & 
Schaeffer 1989, Bernardeau & Schaeffer 1992, Bernardeau & Schaeffer 1999), 



f°° dy r(l + %-<Xy)l t ,« f°° dy . , r(l +5)y- My) 



(31) 



The function 4>(y) (t(jj)) plays an important role in any calculation involving hierarchical ansatz because it completely determines the behavior 
of the PDF p(S) (bias b[S) ) for all values of S. The different asymptotic expressions of <f>{y) govern the behavior of p(S) for different intervals of 
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Figure 3. The skewness parameter S3 and the kurtosis parameter S£ is plotted for different redshift bins. Three different cosmologies are displayed as before. 



5. For large y we can express 4>(y) as: 4>(y) = ay ~ w . No theoretical analysis has been done so far to link the newly introduced parameter to and 
the initial power spectral index n. In the highly nonlinear regime, numerical simulations are generally used to fix lo for a specific initial condition. 
(Colombi et. al. (1992, 1994, 19951, iMunshi et all dl999l) ). Typically for a power law initial power spectrum with pictorial index n = —2 one 
obtains tj = 0.3. The VPF 4>(y) and its two-point analog r(y) both exhibit singularity for small but negative value of y 3 , 

-1 



<t>{y) =4>a- a s T(u s )(y - y s y 



T iv) 



bs(y - y s y 



(32) 



For the factorisable model of the hierarchical clustering the parameter lo s takes the value —3/2 and a 3 and b s can be expressed in terms of the 
nature of the generating function G(r) and its derivatives near the singularity r s (Bernardeau & Schaeffer 1992): 



r(-i/2) 



G'{Ts)Q"(Ts) 



2G'(T s )g"(T s )] 3/2 



Q'"{t 3 ) 



2G'(r s )g"(r 3 



li/2 



Q"'{t s ) 



(33) 



As mentioned before the parameter k a which we have introduced in the definition of Q (r) can be related to the parameters a and to appearing in 
the asymptotic expressions of (j)(y) (Balian & Schaeffer 1989, Bernardeau & Schaeffer 1992), 



u = k a /(k a + 2), 



(34) 



Similarly the parameter y s which describes the behavior of the function 4>{y) near its singularity can be related to the behavior of Q(t) near r s 
which is the solution of the equation (Balian & Schaeffer 1989, Bernardeau & Schaeffer 1992), r s = G' (t s )/Q" (t„), finally we can relate k a to 
y s by following expression (see eq. d34ll): y s — —t s /G'(t s ), or we can write: 



1 (fc a + 2)*°+ 2 

k a (fe a +l) fe » + 1 ' 



(35) 



The newly introduced variable x+ will be useful to define the large S tail of the PDF p(8) and the bias b(S). The asymptotes of (j)(y) are linked with 
the behavior of p(S) for various regimes of 8. For very large values of the variance £2 it is possible to define a scaling function p(5) = h(x)/xi% 
which will encode the scaling behavior of the PDF, where x plays the role of the scaling variable and is defined as x = (1 + 8) We list below 
different ranges of 8 and specify the behavior of p(S) and b(S) in these regimes (Balian & Schaeffer 1989). 



>l + <5>£; p{5) = 



a 1 



1 + 8 



II r(w) v & 



b(8) = 



i» 



LO S 1/2 

2~a) r[i(l + w)] 



1 + 8 



(l-w)/2 



a s /l + 8\ j l + 8\ 

72 [ -~r~ ) ex P ( T ) ' 



6(5) 



(1 + 5) 



..2 ; 



52 



(36) 



(37) 



The integral constraints satisfied by scaling function h(x) are S% = / °° xh(x)dx — 1 and 5j = J °° x 2 h(x)dx = 1. These take care of normal- 
ization of the function p(S). Similarly the normalization constraint over b(S) or equivalently b(x) can be expressed as = xb(x)h(x)dx = 
1, which translates into dSb(S)p(S) = and J x ± dSb(5)p(5) — 1. Several numerical studies have produced the behavior of h(x) and b(x) for 
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Figure 4. The left panel shows the la error in estimation of the variance (k 2 ) c and the right panel shows the error of estimation in the skewness paramter 53 as a 
function of tomographic redshift. The angular scale is fixed at 8 S = 5' . Various curves correspond to a choice of the intrinsic ellipticity distribution of galaxies u f 
and number of galaxies n 3 (=number of galaxies/arc min 2 ) are as depicted. The solid lin es are the difference between the ACDM model and the qCDM model. The 
scatter is computed using the formalism developed in iValag eas. Munshi & Barber 20051). 



different initial conditions (e.g. Colombi et al. 1992,1994,1995; Munshi et al. 1999, Valageas et al. 1999). For very small values of 8 the behavior 
of p(S) is determined by the asymptotic behavior of cj>(y) for large values of y, and it is possible to define another scaling function g(z) which is 
completely determined by ui, the scaling parameter can be expressed as z = (1 + 5)a _1// ^ 1_w ^H' / ' . However numerically it is much easier 
to determine uj from the study of <r(y) compared to the study of g(z) (e.g. Bouchet & Hernquist 1992). 

To summarize, the entire behaviour of the PDF p(8) and bias b(8) are encoded in two different scaling functions, h(x) and g(z). These scaling 
functions are relevant for small and large 8 behavior of the function p(S) and b(S). Typically the PDF p(S) shows a cutoff at both large and small 
values of 8 and it exhibits a power-law in the middle. The power law behavior is prominent in highly non-linear regime. With the decrease in £2 the 
range of 8 for which p(S) shows such a power law behavior decreases finally to vanish for the case of very small variance i.e. in the quasi-linear 
regime. Similarly the bias is a very small and slowly varying function for moderately over dense objects but increases rapidly for over-dense 
objects. These deductions are in qualitative agreement with results from the halo model based approaches. 



5.2 The Quasi-linear Regime 

The Generating function formalism was used by teernardeaul 1 9921 . 1 1 994l) in the quasilinear regime. Unlike highly nonlinear regime the quasilinear 
regime can be dealt with using perturbative analysis. This particular analysis assumes the variance of the smoothed density contrast is smaller than 
unity. The generating function formalism was use to construct the PDF and bias using the tree-level perturbation theory to arbitrary order. In this 
regime the scaling parameters ui and k a can be expressed in terms of the initial power spectral index n. In general the numerical values of the 
parameters k a or uj characterizing VPF or CPDF are different from there highly non-linear values. The PDF and bias now can be expressed in 
terms of Gs(t) (Bernardeau 1992; Bernardeau 1994): 

V ' -G' s (t) L 27T& J V 2£ 2 J ' W 

&(t) = G(t) - 1 = 5. (40) 

The above expression is valid for 8 < 8 C where the 8 C is the value of 8 which cancels the numerator of the pre-factor of the exponential function 
appearing in the above expression. For 8 > 8 C the PDF develops an exponential tail which is related to the presence of singularity in <p(y) in a 
very similar way as in the case of its highly non-linear counterpart (Bernardeau 1992; Bernardeau 1994). 

rx\jx 3a s ^& 5/2 r , <5 \4>s\~\ , x 1 (1 + 5) ,,,, 

p(8)d8= 8 1 exp - \y s ■=- + ^ \d8; b(6) = -— j-^- — = — (41) 

These expressions w ere used bvlMunshi & Jainl d200lL [200ch and Valageas (2000) for the construction of weak lensing PDF in projection. The 
bias was studied in dMunshil 120001) for projected surveys. The tests against numerical simulations show remarkable agreement for a range of 



[(i + foW) 



1/fea 



(39) 
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Figure 5. The PDF p(»?) of the reduced convergence r\ as a function r\ = 1 + 8. The plots with decreasing peak height correspond to lower redshift bins. Two different 
approximations are being compared. The solid line correspond to the lognormal approximation and the dashed line correspond to the perturbative calculations. The 
results are shown for a smoothing angular scale 8 S = 5'. 



angular scales. Later studies refine d these results as well as incorporated various different smoothing windows I Munsh i. Valageas & Barberil2004 
I Valageas. Munshi & B arber 2005). We extend these results derived for surveys in projection to tomographic surveys in this paper. 

It is worth mentioning that, there have b een various attempts to extend the perturbative results to the highly nonlinear regime (see e.g. 
JColombi et all 19971 : Ivalageas & Munshill2004ri ) 



6 THE LOGNORMAL DISTRIBUTION 



An alternative to the hierarchical ansatz, the lognormal distribution, for the description of the matter PDF which has long been known as a suc- 
cessful empirical prescription for th e characterization of the dark matter distribution as well as the observed galaxy distribution lHam iltoi][l985r . 
IColes & Jones|[l99ll :lBouchet et al 1993; Kofman et al. 1994). Detailed discussion for comparison of lognormal distribution and the perturbative 
calcu lations can be found in teernardeau & Kofmanlll 99*51) . The lognormal distribution was further generalized to the sfcewerf-lognormal distribu- 
tion jColombill99"4l) . In general a variable might be modelled as lognormal if it can be thought of as the multiplicative product of many independent 
random variables. 

Although inherently local in nature, the lognormal distr ibution can provide a good fit to both one-point PDF and its gener alisation to compute 
its tw o-point analog and hence the bias jTaruva et al j|2002h . The one- and two-point lognormal PDF can be expressed as jKavo. Taruva. Sutol 
1200 lh : 



p ln (S)d8 



/2ttE 

Pln(Sl,S2)d6ldS2 = 



: exp 



"2S 2 
1 



X 12 



d8 
1 + 5' 

exp 



ln(l + a 2 ); 



E(A 2 



2X12A1A2 



ln[(l + 5)^(1 + ff 2 ); 
dSi d&2 



2(E 2 



1 + Si 1 + S 2 ' 



Aj = ln[(l + 8i)\J (1 + II); Ei 2 = ln(l + 62) 
In the limiting case of large separation X12 — > we can write down the two point PDF 

Pln(8l,8 2 ) = Pln(<5l)pin(<52)[l + 6ln {Si )£l2&ln {S 2 )] J 



bin($i) = Ai/E». 



(42) 
(43) 
(44) 

(45) 



It is however easier to estimate the cumulative or integrated bias associated with objects beyond a certain density threshold <5o- This is defined 
as bin{8 > 80) = p\n{S)bi n (S)dS / p\ n (8)dS. In the low variance limit £2 — > the usual Gaussian result is restored b(S) = 5/^2. The 
parameters A, Aj, X12, E that we have introduced above can be expressed in terms of the two-point (non-linear) correlation function £12 = {8182) 
and the nonlinear variance a 2 = (8 2 ) of the smoothed density field. To understand the construction of lognormal distribution, we introduce a 
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Gaussian PDF in variable x\ p(x) — (2tyT, 2 ) 1 ^ 2 exp[— (x — p) 2 /2E 2 ]. With a change of variable x = ln(t) we can write down the PDF of y 
which is a lognormal distribution p(t) = (27rE 2 ) 1,/2 exp[— (ln(t) — p) 2 /2Y?]/t. The extra factor of (1/x) stems from the fact: dt/t = dx. Note 
that y is positive definite and is often associate with p/po = 1 + 5 which means (t) = 1. The moment generating function for the lognormal in 
terms of the mean p and the variance E has the following form: (t n ) = exp(np + n 2 E 2 /2). This however leads to the fact that if the underlying 
distribution of x or the density is Gaussian we will have to impose the condition: p = — E 2 /2. Here in our notation above E is the variance of 
the underlying Gaussian field. The variance of t defined as (t 2 ) — (t) 2 = exp(E 2 ) — 1 = a 2 . So we can write E 2 = ln(l + a 2 ). This is the 
result that was used above . The generalization to two-point or bi-variate PDF can be achieved following the same arguments and can be found in 
dKavo. Taruva. Sutdl200ll) . 

In contrast to the lognormal model the widely-used non-local ansatz for the evolution of the variance that was introduced by dHamilton et al.l 

1985), the evolved nonlinear two-point correlation function is linked to that of the initial or linear two-point correlation at a different scale: 
?2 ! (-Rm) = fnl{$2 n (Rim)]- The two different length scales are related by the following expression: = (1 + |j ! {Rm))Rfi n - Such as an 
ansatz is derived using pair conservation equation. The non-linear (Eulerian) length scale R n \ is linked to the linear (Lagrangia n) length scale Rn n 
from where the structure has collapsed. Numerical simulations are typically used for the determination of th e fitting function / n i {P eacock & Doddj 
1 19941) . However, simpler asymptotic power-law forms exist in different regimes of gravitational clustering dMunshi &^ Padmanab hanll 19971) . 

The validi t y and limita tions of the one-point and two-point PDFs have been studied extensively in the literature against N-body simulations. 
In lBernardeav] dl992l . 1 19941) it was shown that the PDF computed from the perturbation theory in a weakly nonlinear regime approaches the 
lognormal distribution function only when the primordial power spectrum is locally of the form P(k) oc k n " with the effective local spectra slope 
of the power spectrum n e ~ —1. It was also shown that in the weakly nonlinear regime the lognormal distribution is equivalent to the hierarchical 
model with a generating function G(t) = exp(~r). This leads to the following skewness and kurtosis parameters: 

S% = 3 + a 2 ; S2 = 16 + 15cr 2 + 6<r 4 + cr 6 . (46) 

In general the a 2 — > leads to Sp = p p ~ 2 - On this basis iBernardeau & Kofmanl dl995h argues that the agreement of lognormal PDF with 
numerical simulations should be interpreted as purely accidental and the success of the lognormal model is simply related to the fact that for all 
scales relevant to cosmology the CDM power spectrum can be approximated with a power law with effective slope n e ~ — 1. However subsequent 
studies using numerical simulation it was shown by various authors that the lognormal distribution very accurately describ es the cosmological 
distrib ution functions even in the nonlinear regime a < 4 for a relatively high values of density contrast 5 < 100 (see e.g. iKavo. Taruva. Sutq 
J200lh ). 

There is no complete analytical description of gravitational clustering in the highly nonlinear regime. However several dynamical approx- 
imations were proposed in the past to mimic certain fe atures of gravitational clustering neyond the weakly nonlinear clsutering. The Frozen 
Flow App roximation (FFA) approximation proposed by dMatarresse et alJI 19921) is one such approximation, using perturbative techniques it was 
shown by dMunshi. Sahni. Starobinsld[T994l) that the FFA develops exactly the same generating function as the lognormal approximation in the 
quasilinear regime. 

The error estimates for various lower o rder S v (right panel) and the varia nce (left-panel) are shown in Figure {4]l. A complete analytical 
formalism for calculation of error are given in Valageas, Munshi & Barbed d2005h . It requires the knowledge of higher order S p parameters. In our 
calculation of error the higher order S p parameters are modelled according to lognormal distribution. Different levels of noise as chareterized by 
the parameters that describe intrinsic ellipticity distribution e and number density of galaxies/arcmin 2 or n g are considered. 



7 THE PDF AND BIAS OF SMOOTHED REDSHIFT-RESOLVED CONVERGENCE MAPS 



For computing the probability distribution function of the smoothed convergence field for individual tomographic maps k w (9o), we will begin by 
constructing its associated cumulant generating function for individual tomographic bins 3>i+ K (g n (y)- The construction is based on modelling of 
the volume-averaged higher order correlation function (k^ (#o))c in terms of the matter correlation hierarchy: 

(y)=v + T ( *W (g ° ))c V P (47) 

p=2 \ K (i)\ a °))c 

Now using the expressions for the higher moments of the convergence k(8o) in terms of the matter power spectrum, Eq[l0]and Eq|20|gives: 

SiWW-l. pi y -y + J Z. p! S P dA(r)2 ( P -i) (r) [^) K(eo) J V' $2 -<*(*o)>c. (48) 
We can now use the definition of cj>(y) for the matter cumulants to express &i+ K (o ) (y), m terms of (j)(y): 



, \ r A r ^gW lJ gwW tool?) 1 min min f , , s 
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Figure 6. The cumulative bias b(> r]) of the reduced convergence r) = (1 + S) is plotted as a function r\ for various redshift bins. The smoothing angular scale 
is 9q = 5'. As before two different approximations are considered. The lognormal approximation (solid lines) and the perturbative calculations (dashed lines) 
reproduce nearly identical results. The curves that saturates at a higher values of cumulative bias for higher values of rj correspond to larger smoothing angular 
scales. 



The extra term comes from the p = 1 term in the expansion of $ 1+fc (e ) . Note that we have used the fully non-linear generating function qS for 
the cumulants, though we will use it to construct a generating function in the quasi-linear regime. The analysis becomes much easier if we define 
a new reduced convergence field: 

Vd)(Qo) = (n(l] n - «;(i)(6>o))/K(?) n = 1 + K(i)(0o)/|K(i) m |- (50) 

Here the minimum value of K(i) (6o) i.e. t™j n occurs when the line-of-sight goes through regions that are completely empty of matter (i.e. 5 = — 1 
all along the line of sight in a redshift window that defines a specific bin i.e. Eq.([3j). While Ku) (6q) depends on the smoothing angle, its minimum 
value «™j n depends only on the source redshift and background geometry of the universe and is independent of the smoothing radius. With the 
reduced convergence r\, the cumulant generating function is given by, 

lK! [KfflJa Je (r) rld\(r) g)-( do) y \ 

The thus constructed cumulant generating function & v (y) satisfies the normalization constraints SJ 7 = = 1. The scaling function associated 
with P(rf) can now be easily related with the matter scaling function h(x) introduced earlier: 

ft« (a0 = _f ^L exp{xy) ^ {y) . h <S) {x) = 1 r dr \ J^pl. ]\{ . Y (52) 

While the expressions derived above are exact, and are derived for the most general case, using only the small angle approximation, they can be 
simplified considerably using further approximations. In the following we will assume that the contribution to the r integrals can be replaced by 
an average value coming from the maximum of U)u)(r), i.e. r c (0 < r c < r 3 ). So we replace J f(r)dr by l/2/(r c )A r where A r is the interval 
of integration, and f(r) is the function of comoving radial distance r under consideration. Similarly we replace the w(r) dependence in the k 
integrals by oju) (r c ). 



I mm, _ 



^W (rc) ' K 1 ^ [ J J^P^)WMm)\ . (53) 

Under these approximations we can write: (y) = (f>^ (y) and (x) — h(x) . Thus we find that the statistics of the underlying field rj = 1 + 5 
and the statistics of the reduced convergence r\ are exactly the same under this approximation. Though we derived the results from considering 
a specific form of hierarchical ansatz, the final result is remarkably general. It simply means that independent of detailed modelling the reduced 
convergence will always follow the statistics of underlying mass density distribution. 

Finally, we can express the relations connecting the probability distribution function for the smoothed convergence statistics k(6o), the 
reduced convergence r)(9o), i.e. for individual bins we can write: 

p (i) («)=p(»7 Ci) )/|«« in |. (54) 
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Figure 7. The left panel shows the PDF of the redshift resolved convergence and the right panel shows the associated cumulative bias. The smoothing angular scales 
considered is 9q =5'. Only three tomographic bins are chosen for display to avoid cluttering. Two different approximations are used; the lognormal distribution and 
the perturbative calculations. The approximations give near identical results. Three different redshift bins are displayed z s = 0.698, f .095, f .493 (top, middle and 
bottom curve respectively). 



This is one of the most important results in this paper. A few comments are in order, it is possible to integrate the exact expressions of the scaling 
functions, there is some uncertainty involved in the actual determination of these functions and associated parameters such as aj,fc a ,:r* from 
N-body simulations and its unclear how much is there to gain by doing exact calculations that involve approximate picture for the underlying mass 
distribution; see e.g. Munshi et al. (1999), Valageas et al. (1999) and Colombi et al. (1996) for a detailed description of the effect of the finite 
volume correction involved in their estimation. Throughout our an alysis we have used a top-hat fi lter for smoothing the convergence field, but the 
results can easily be extended to compensated or Gaussian filters dBemardeau & Val ageas 2000). Using Eq ,d54t one can derive an approximate 
expression for the lower order moments which can, in turn, be used to derive order of mag nitude relations S* = S^/[k^ n ] p - 2 for p > 2. The 
parameters Sp for the underlying density contrast 8 are specified by t he choice of a specific hi erarchical model. The computation of error bars for 
these lower order moments can be done using formalism developed in lMunshi&Colesl l l2003l) . 

In Figure (|5j we show the PDF of the reduced convergence r\ for smoothing angular scale 9o = 5' for the lognormal and hierarchical 
approximation discussed above, for individual redshift bins as well as for a projected survey. 



7.1 The bias associated with convergence maps 

To compute the bias associated with the peaks in the smoothed convergence k field we have to first develop an analytic expression for the generating 
field /3 K (j/i, y-2) for the convergence field k(0o)- We will avoid displaying the smoothing angular scale 6*o for brevity. Throughout the statistics are 
for smoothed convergence fields. For that we will use the usual definition for the two-point cumulant correlator C pq for the convergence field (for 
a complete treatment see Munshi & Coles, 1999b): 

= (^(ntf^^^^c/KTr 1 ^] 9 - 1 ^. (55) 

We will show that, like its density field counterpart the two-point generating function for the convergence field k can be expressed (under certain 
simplifying assumptions) as a product of two one-point generating functions /3 (y) which can then be directly related to the bias associated with 
"hot-spots"in the convergence field. 

ti ij) (yi,v*) = E = E ^E^' = ^(yM j) (v2). (56) 

p,q P 9 

It is clear that the factorization of generating function depends on the factorization property of the cumulant correlators i.e. Cp q = C^C^. Note 
that such a factorization is possible when the correlation of two patches in the directions Cli and 0,2 [£12 ^] is smaller compared to the variance 
[£vw] for the smoothed patches We will now use the integral expression for cumulant correlators (Munshi & Coles 1999a) to express the generating 



© 0000 RAS, MNRAS 000, 000-000 



14 Munshi et al. 



ACDM vs Dark Energy Models 



20 



10 



-10 



^~ 




wOwl— A 



: U \ 




\\ o =5- : 






wOwl-A V N 's. '_ 




, i i 1 ,\, \ 



-0.02 0.02 0.04 



-0.02 



0.02 0.04 



Figure 8. We compare the PDF and cumulative bias associated with tomographic convergence maps for various tomographic bins for two different dark energy 
models against the ACDM model. The left panel plots the ratio A p (k) = p(re) — Pacdm(k) for a smoothing angular scale of 6o = 5', The right panel depicts 
the ratio (> ft) = b«(> ft) - 

'acdm^ k ) ^ or ^ cumulative bias. Three different redshift bins are displayed z s = 0.698, 1.095. For a given smoothing 
angular scale and a fixed redshift, the (thick) curves with higher positive peak heights correspond to the model too = —1,10 i = 1 and the (thin) ones with lower 
peak heights correspond to uiq = —0.9, w\ = 0. 



function which in turn uses the hierarchical ansatz and the far field approximation as explained above 

°° fv(ij) 1 1 i rr, (jjP (r)u> q (r) 

^^^ikw^m^mL dr 4,)^-') [ ^- (0][ ^ (r)1P+9 " 2 ^- (57) 

It is possible to further simplify the above expression by separating the summation over dummy variables p and q, which will be useful to establish 
the factorization property of the two-point generating function for bias ,1/2). We can now decompose the double sum over the two indices 

into two separate sums over individual indices. Finally, using the definition of the one-point generating function for the cumulant correlators we 
can write: 

The above expression is quite general within the small approximation and large separation approximations, and is valid for any given specific 
model for the generating function Qsij). However it is easy to notice that the projection effects as encoded in the line of sight integration do not 
allow us to write down the two-point generating function /3 K (yi, 3/2) simply as a product of two one-point generating functions /3 v (y) as was the 
case for the density field 1 + 5. 

K ' Jo [«ft m ] K5) n l [& )K ] M^ ,)K ] d iM l K (i) in l Mct )K ] <W l«w n V 

We use the follwoing equation in assocition with Eq.J53b to simplify the above expression: 

Use of these approximations gives us the leading order contributions to these integrals and we can check that to this order we recover the fac- 
torization property of the generating function i.e. (3/1, 2/2) = Pri\yi) Pt)\v2) = /? 1 ^(yi)/3 1 +j(j/2). So it is clear that at this level of 
approximation, due to the factorization property of the cumulant correlators, the bias function b v (x) associated with the peaks in the convergence 
field k, beyond certain threshold, possesses a similar factorization property too as its density field counterpart. Earlier studies have established 
such a correspondence between convergence field and density field in the case of one-point probability distribution function p(S) (Munshi & Jain 
1999b), 

6«(xi)ft«(xi)6«( :C2 )ftV' ) (x2) = b^ t (xi)h^l f (x 1 )l^l s (x2)h^. s ( i X2), (60) 
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where we have used the following relation between /3 v (y) and b v (x). For all practical purpose we found that the differential bias as defined above 
is more difficult to measure from numerical simulations as compared to its integral counterpart where we concentrate on the bias associated with 
peaks above certain threshold. The cumulative bias b v (> x) can also be defined in an analogus manner: 

b v (x)hr,(x) = -— dyr(y)exp(xy); b n (> x)h v (> x) = -— / dy exp(xy). (61) 

2" J_ ioo 2m J_ iaa y 

It is important to notice that, although the bias b(x) associated with the convergence field k and the underlying density field are identical, the 
variance associated with the density field is very high, while projection effects substantially reduce the variance in the convergence field. This 
indicates that we have to use the integral definition of bias to recover it from its generating function; see Ea.d61b- Now writing down the full two 
point probability distribution function for two correlated spots in terms of the convergence field k(9q) and its reduced version r\: 

P ^^K 1 ,K 2 )dK 1 dK 2 =P^\K 1 )p^\ K2 )(l + b (i \ Kl )^b U) (K2))dK 1 d K2 , (62) 

Pv(Vi>V2) dm dr)2 = P ( v l) (Vi )Pv ) (V2 ) ( 1 + b^ ( m )tf 2 b v J) (V2)) dm drj 2 . (63) 
Using Eq.d54t that p^ (k) = p v {ri) /l K mLl we a ^ so not i ce mat £12^ ~ £i 2 /[ K mL] [ K min]' fr° m which we can now write: 

6 W (k)=&«(»/)/|k£LJ. ( 64 ) 

Together with Eq.J54t. Eq.ll64t can be used to construct analytical estimates of pdf and bias about individual bins. In addition these results are 
applicable to the modelling of joint PDFs involving two separate redshift bins. 

Figure (|6) shows the bias associated with the reduced convergence for individual bins as well as for the entire survey. The smoothing angular 
scale is 9q = 5'. In Figure (0 shows the PDF and bias associated with the convergence k. In Figure ([8} we plot the difference in PDF between 
various cosmological scenarios. 



8 EFFECT OF NOISE ON ONE- AND TWO-POINT PDF 

The PDF we have considered so far are free from noise. In this section we will present the resu lts of estimates of error relating to the PDFs, those 
at the level of one-point and two-point PDFs. T hese results will generalise the ones found bv lMunshi & Colesl d2003l) for lower order moments 
and later by Valageas, Munshi & Barbed d2005h . Inclusion of noise can be incorporated through a convolution. We will assume the noise to be 
Gaussian and uncorrelated with the signal. However the variance of the noise will depend on the surface density of galaxies in individual bins. 
With these simplifying assumption, for the i-th tomographic bin we can write: 

/oo 
pW(K-n)p$(n)dn. (65) 
- 00 

Here Pq (n) is the noise PDF assumed Gaussian, and pn ( K ) is the convergence PDF in the presence of noise (the subscript G denotes Gaussian). 
We take — of /(2n s 7r#o)- Here cr E is intrinsic ellipticity distribution of galaxies and n 9 is the number denisty of galaxies and 60 is the 
smoothing angular scale. The two-point PDF can also be modified to include the effect of noise in a similar manner. The equivalent expression for 
2PDF can be written as: 

p^' ) (« 1 ,« 2 )=pl i) (K 1 )p^ ) («2)(l + 6l i H«i)^2#(«2)) J (66) 
Which is obtained by convolving the noise PDF with the 2PDF: 

/oo 
- n,i,«2 - n 2 ) Pq (ni) p% {n 2 ) dm dn 2 . (67) 
- OO 

Comparing Eq.d66t and Eq.J67li we can write the expression for the noisy bias b n {n) as: 

/oo p 00 

p {i) (k - n)fo (l) (k - n)p { § (n)dn / I p {l) (k - n)p { § (n) dn. (68) 
- 00 J — 00 

Notice that depending on redshift distribution of sources, the noise maps n^' and n^' can be different for two tomographic bins. We also assumed 
that noise in different bins are statistically independent. The cumulative bias for the i-th tomographic bin that include noise b„\> k) can be 
expressed in terms of pn («) («) just as its noise-free counterpart: 

poo roc 

^ ) (>«)=/ 9$(K)b®{K)dk/ I pttH^dk. (69) 
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Figure 9. The sacatter in estimation of binned PDF for a given intrinsic ellipticity distribution and sky coverage is displayed. We assume an all-sky coverage. The 
effect of intrinsic ellipticity distribution is displayed in the left panel. The right panel depicts the effect of binning. The angular scale in each case is 9o = 5' and the 
redshift is z s = 0.698. The expression for <t(p(k)) is given in Eq.{72j. For the left panel we consider A = 0.02 and for the right panel n g = 25. A background 
ACDM cosmology is accumed for these calculations. 



Errors associated with binned tomographic noisy PDF can be analysed using following quantites: 



N = n„7r<95 = 314 9 - (70) 

V 100 arcmin / \ 1 arcmm J 

Here n g is the number density of galaxies, 9 3 is the smoothing angular scale in arc-minutes for a given survey strategy. For a given survey we 
denote the area covered by A and introduce a parameter N c which will be used in expressing the signal-to-noise estimates of the PDF p(tt). We 
define the following variable that will be useful in qunatifying scatter in a noisy PDF: 

Nc = ^_ = 2 . 7 x 10 *C d )(^L_Y\ (71) 

(26»o ) 2 V. 300 degree 2 / V 1 arcmm ^ 



Finally the scatter <j(p(/t)) in the measured convergence power sepctra p(k) can be expressed as (Vala geas. Munshi & Barbe3l2005h : 



\ 1/2 

l 



- 1 



(72) 



These expression can be modified and used to compute the scatter in individual redshift bins by simply changing n g to surface density of individual 
bins fig and p(k) to p' 1 '(k). The source density for individual bins for a given survey can be computed using = J^ z+1 p s {z)dz. The bin 

width A is left as a free parameter. In Figure l[9} we have plotted the scatter a(p„ («) as a function of intrinsic ellipticity distribution and bin 
width for a smoothing angular scale of 9o = 5' and redshift z s = 0.698. The results for the difference in noisy PDFs are plotted in Figure l |10t . 



9 CONCLUSIONS 

Previous tomographic studies of weak lensing have typically worked with the lower order cumulants; we have generalized here these results to the 
case of the entire one- and two-point PDF, which contain information about the cumulants to an arbitrary order. The construction was performed 
using a generating function formalism based on hierarchical ansatz and a lognormal model. Our analysis generalizes previously obtained results 
derived for projected or 2D surveys. Though we have considered a top-hat filer c onvergence maps due to their simplicity, similar results can be 
obtained for related statistics such as the shear components or aperture mass M ap dBemardeau & V alageas 2 000h . 

The PDFs for the individual bins are constructed by generalization of the previously introduced global variable K mm , for individual bins 
i.e.K mm , that was used in the context of 2D projected maps. Next, using K™j n , reduced variable t]^' is defined for each individual bins whose 
statistics can directly be linked to that of underlying density contrast 8. The convergence in individual bins can then be mapped to unique values 
of r] = 1 + 8 for a given smoothing angular scales 6*o. 

For modelling the statistics of underlying density contrast 8 we have assumed two completely different model: the hierarchical ansatz along 
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Figure 10. The difference of noisy ACDM PDF and dark energy models A p (re) = p(re) — PACDM (re) is plotted as a function of re. The smoothing angular scale, bin 
size and galaxy number density is as depicted. The scatter in estimation is smaller compared to the difference in the PDFs considered. The comsological parameters 
considered are the same as the ones in Figure-{8). The two survey configurations that we have considered both produces near identical results. 



with its perturbative counterpart as well as the lognormal distribution. Both these approximations have been used successfully in various cosmo- 
logical contexts. There are a wide class of models that are available under the general category of hierarchical ansatz. The main motivation behind 
our choice of a specific hierarchy is simplicity. In recent years more sophisticated models of hierarchical clustering have been proposed which 
rely more on inputs from numerical simulations. The generic results we have derived here indeed ca n be improved us i ng su ch modelling though 
the fundamental characteristics will remain unchanged. In our treatment we find, in agreement with lMunshi & Wand d2003h . the dynamical and 
geometrical contribution can be treated separately. The geometrical effects are completely encoded in a parameter n mm . The reduced convergence 
as defined is independent of the background geometry of the universe and essentially probe the evolution of gravitational clustering. We showed 
that a set of ^™) ln defined for a given set of redshift slices are adequate to characterize not only individual PDFs for each bin but it is also sufficient 
to study the joint two-point PDF among two different bins. The PDF of the reduced convergence rf % ' for individual bins or joint PDFs for a pair 
of bins generalizes the earlier studies where the projected or 2D maps were considered in a straight forward manner. 

We also note that the construction of convergence maps is difficult compared to the direct evaluation of non-Gaussian statistics from shear 
maps. On the other hand convergence statistics can directly be modelled at arbitrary order whereas for shear field the computation is done mostly 
order by order manner. An independent analysis of convergence maps constructed from shear maps should therefore be useful in constraining 
various errors which might get introduced during various stages of data reduction. 

In our analysis we have ignored the noise from intrinsic ellipticity distribution of galaxies as well as from shot noise res ulting from finite 
numbe r of galaxies that are used to t race t he underlying mass distribution. These issues have been dealt with in great detail in lMunshi & Colesl 
J2003b : Ivalageas. Barber & MunshJ fe004) . Dividing the source population into bins reduced the number-density of sources. This in turn will 
increase the level of noise or the scatter in the estimator. In our analysis we have considered two different survey configurations, i.e. LSST and 
DES and found that for our choice of tomographic bins the one- and two-point PDFs are very similar in nature. 

The lognormal distribution has already been use d to model the statistic s of weak lensing observables jMun shi 2000; T aruva et alj|2002l) and 
the clustering of Lyman alpha absorption systems e.g. tei&Davidsonl 19971) . One-to-one mapping of initial density fields to evolved density fields 
using maps that are consistent with lognormal distribution function was not found to be very successful and the success of a lognormal distribution 
function in reproducing the statistics of gravitational clustering still remains somewhat unclear. 

Tomographic weak lensing surveys can be cross-correlated with external data sets including frequency cleaned maps of secondaries from 
ongoing CMB surveys; e.g. the thermal Sunyaev-Zeldovic (tSZ) maps or y-maps that will be available from CMB surveys such as Planck. The 
cross-correlation with tomographic information can help to understand the evolution of cosmological pressure fluctuations responsible for tSZ 
effect with redshift. The formalism presented here is perfectly suitable for such analysis. Detailed results of such analysis will be presented 

elsewhere. In addition to the weak lensing surveys the Supernova pencil beam surveys might also benefit for the results presented here. 

To summarize, we have exte nded results derived in three different previous papers ( Munshi & Colej2003l : IValageas, Munshi & Barbeill2005l ; 

Valageas, Barber & Munshi 2004) to take into account tomographic bins within which the photometric redshift are available. The results obtained 
previously for one-point PDF are now extended to two-point PDF. These results can provide an alternative to usual Fisher-matrix analysis that 
is employed to optimize survey strategies. We have concentrated mainly on analytical results in this paper. The numerical results regarding 
optimization of survey strategy using these results will be considered elsewhere. 
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